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ABSTRACT 

We present a novel, virtual X-ray observatory designed to obtain synthetic observa- 
tions from hydro-numerical simulations, named PHOX. In particular, we provide a 
description of the code constituting the photon simulator and of the new approach im- 
plemented. We apply PHOX to simulated galaxy clusters in order to demonstrate its 
capabilities. In fact, X-ray observations of clusters of galaxies continue to provide us 
with an increasingly detailed picture of their structure and of the underlying physical 
phenomena governing the gaseous component, which dominates their baryonic content. 
Therefore, it is fundamental to find the most direct and faithful way to compare such 
observational data with hydrodynamical simulations of cluster-like objects, which can 
currently include various complex physical processes. Here, we present and analyse 
synthetic Suzaku observations of two cluster-size haloes obtained by processing with 
PHOX the hydrodynamical simulation of the large-scale, filament -like region in which 
they reside. Taking advantage of the simulated data, we test the results inferred from 
the X-ray analysis of the mock observations against the underlying, known solution. 
Remarkably, we are able to recover the theoretical temperature distribution of the two 
haloes by means of the multi-temperature fitting of the synthetic spectra. Moreover, 
the shapes of the reconstructed distributions allow us to trace the different thermal 
structure that distinguishes the dynamical state of the two haloes. 
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1 INTRODUCTION 

As optimal laboratories for both cosmology and astro- 
physics, galaxy clusters have been thoroughly investigated 
with large surveys as well as via dedicated single-object ob- 
servations. In this respect, X-ray measurements still rep- 
resent one of the best ways to fully study their struc- 
ture, traced by the hot plasma filling their potential wells. 
In particular, the most important and difficult quantity 
to infer is the total gravitating mass, which can be de- 
termined from measurements of the intra-cluster medium 
(ICM) temperature and density profiles, assuming hydro- 
static equilibrium to hold. So far, studies of the X-ray- 
emitting ICM have been mainly limited to the innermost 
region, where the investigation of the baryonic physics is 
observationally less challenging than in the outskirts. Al- 
though very difficult to map in the X rays, a character- 
ization of the outer region of clusters is in fact crucial 
to understand the formation and evolution of this struc- 
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tures and to use them as cosmological probes, since it en- 
closes a significant fraction of the cluster volume, where 
the properties of the accreting ga s and of the dark matte r 
halo are still uncertain (e.g., see lEttori fc Molendil lioilh . 
Only very recently, X-ray observations performed with the 
Japanese satellite Suzaku ha ve reached the virial radius for a 
few clusters of galaxies (e.g. Fuiita et al.ll2008l ; iBautz et al 



2009| : lGeorge et alj2009l;lReiprich et al.ll2009l;lHoshino et al 



20101: iKawaharada et "all |2010t iSimionescu et all 120111 : 
Akamatsu et alj 201 ll ). promisingly taking a step forward on 
the temperature profile debate and indirectly opening new 
perspectives in controlling the biases on cluster-based cos- 
mological investigations. These results are just one example 
of the improvements by which X-ray observations continue 
to provide us with an increasingly detailed picture of galaxy 
clusters, for which a clear interpretation of the underlying 
physical processes is very challenging. 

On the theoretical side, many interesting issues on 
cluster formation and evolution are addressed by means of 
hydro-numerical simulations, which represent a powerful 
tool to investigate in detail the nature of these compli- 
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cated astrophysical objects (e.g. iBorgani fc Kravtsovll2009l . 
for a recent review). Simulations of large cosmological 
boxes that include several massive cluster-like haloes can 
now be performed incorporating not only the dominat- 
ing dark matter component but also baryonic matter, 
in the form of gas- and star-like particles (or cells), 
governed by hydrodynamical processes. With particular 
concern for galaxy clusters, current simulations are not 
only able to account for basic gas hydrodynamics, but 
implement also more complicated models for st a r formation 
from multi-phase medium (e.g. iKatj Il992l; iKatz et al.l 
1 19961 ; ISpringel fc Hernquistl |2003l ; iMarri fc White! l2003h 
and thermal or kinetic feedback from supernovae- 



driven winds (jNavarro fc White 1 19931 ; IScannapieco et ail 



ment, metal and molecule coolin 


g (e.g. 


Mosconi et al.l 


2001; Yoshidaet al.l 120031; iTornatore et al 


12004 2007; 


Scannapieco et al. 2005; iMaio et all 120071. |2010|). thermal 



conduction dCleary fc Monaghan 1999; lJubelgas ct al. 2004: 



Ruszkowski et al. 
20051: 



20ld) AGN fe edback (e.g. ISpringel et al l 



Pi Matteo et all 120051: ISiiacki fc Springel 



2006; 



Siiacki et al. 2007, 2008; Puchwein etail 2008; F abian et ail 



201(1 iDubois et al.l l2010l; iTevssier et all l201ll). c osmic 
rays (|Pfrommer et al. 120071; lJubelgas et al l 120081) and 
magn et ic fields dPhillips fc Mona ghan 1985; Dola g~et al.l 



19991; iBriiggen et all 120051 ; [Pricefc Monaghanl 120051 ; 
Dolag fc Stasvszvnl |2009| ), to name the most important 



effects. 

The ideal achievement would be the combination of 
both these sources of information, by directly comparing 
simulated clusters to X-ray observations of real objects. De- 
voted to this goa l, sophisticated numerical codes, such as X- 
MAS/X-MAS2 llGardini et al.ll2004l; iRasia et al.ll2008l) and 
more recently, XIM (|Heinz fc Bruggenl 120091 ; iHeinz et al.l 
l2010l ). have been developed in the last years in order to turn 
hydro-simulation outputs into mock images with a given X- 
ray telescope. In particular, the X-MAS virtual telescope is 
explicitly designed to process outputs obtained from SPH 
codes like GADGET, whereas XIM is particularly dedicated 
to grid-based hydrodynamical simulations, e.g. pe rformed 
with the AMR code FLASH (|Frvxell et al.l |2000| ). which 
makes them fundamentally complementary. Similarly, both 
X-MAS and XIM use a plasma thermal emission code by 
which they calculate the projected emission associated to 
the gas component in the simulation. The computational 
effort, required to calculate this, is usually reduced by in- 
terpolating among the model spectra externally stored in a 
library of templates. Essentially, the emissivity integrated 
along the chosen line of sight is calculated in terms of flux, 
depending on the properties of all the gas at a given pro- 
jected position. Lastly, the convolution with instrumental 
response is performed and the final photon event file is gen- 
erated. Such virtual telescopes represent an important step 
forward with respect to visualization tools used to produce 
surface brightness maps from simulation outputs, since they 
properly calculate the emission associated to the gas com- 
ponent accounting for the structure of the simulated source 
along the line of sight. Moreover, they are able to include ac- 
curate calculations required to obtain realistic mock X-ray 
observations, for instance convolving the simulated spectra 
with any given instrument response and telescope PSF. In 
fact, it is vital that the synthetic observations produced by 



such simulators are as much as possible similar to the stan- 
dards commonly used by observers, so that the comparison 
can be the most faithful. As such, this advanced approach 
easily leads to a challenging increase of the parameters to 
handle in order to realistically model the plasma emission 
and the detailed three-dimensional structure of the simula- 
tions. 



Analogous in the scope, the virtual telescope presented 
in this paper, PHOX, is also dedicated to convert hydro- 
numerical outputs into mock X-ray observations. However, 
the novelty of our simulator lies in the method adopted, by 
which the spectral emission calculated singularly for each 
gas element in the simulation is immediately converted into 
a discrete sample of photons, collected and stored before 
projecting along any line of sight and convolving with any 
desired instrument. With this strategy we are able to signif- 
icantly gain in computational effort, since it requires to pro- 
cess the original simulation only once, independently of the 
specific synthetic observation to be performed afterwards. 
Furthermore, the guarantee for high spatial and spectral res- 
olution, preserved without dramatically increasing the com- 
putational cost, offers the possibility to anticipate the obser- 
vational achievements of upcoming X-ray missions, such as 
IXO/ATHENA, in which the high-resolution spectroscopy 
will allow us to explore the intrinsic structure of galaxy clus- 
ters through the study of their spectral features. In fact, 
it has been suggested in the last decade that the diagnos- 
tics of the broadening of heavy-ion emission lines in highly- 
resolved X-ray spectra of galaxy clusters could actually pro- 
vide valuable information about the u nderlying structur e 
of the ICM and its velocity field (see IPawl et al.l l|2005h . 
for a preliminary, interesting work in this field), useful for 
the dynamical classification of these objects. The expecta- 
tions for such line diagnostics are related in particular to 
the most prominent emission line in X-ray spectra, namely 
the ~ 6.7keV line from helium-like iron. In fact, the large 
atomic mass of the FeXXV ion significantly reduces the ther- 
mal line broadening, so that the line width turns out to be 
definitely more sensitive to turbulent or bulk g as motions 
l|lnogamov fc Sunvaevl2003l ; ISunvaev et al.ll2003l ). This kind 
of studies can eventually help constraining the detectabil- 
ity of non-th ermal motions that are likely to establish in 
the ICM (e.g. iRebusco et alj|2008 | ) and can compromise X— 
ray mass measurements (IRasia et al. 2006; Fang et aill2009l ; 
lLau et al.ll2009l ; iBiffi et al.ll201ll ). 



The paper is structured as follows: Section [2] is devoted 
to the description of the method implemented in the photon 
simulator and the fundamental units of the code. In Sec- 
tion [3] we apply PHOX to the hydrodynamical simulation 
of a filament -like structure, describing in detail each phase, 
from the photon generation till the mock Suzaku observation 
of two massive galaxy clusters residing in the filament and 
the spectral analysis of the synthetic spectra. Results about 
the recovering of the ICM emission measure distribution of 
the two cluster-like haloes from the multi-temperature fit- 
ting are presented in Sectionf!] Finally, we discuss the results 
obtained, summarise and draw our conclusions in Section [S] 
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2 PHOX: THE TECHNIQUE 

The novelty of our photon simulator arises from the method 
implemented. In this Section we describe the three main 
parts of the code that constitute the virtual telescope, 

phosQ 

for which a schematic overview is also presented 
in Fig. 12. II PHOX consists of a C code, which m akes use 
of th e publicly available X-ray package XSPEC (see I Arnaudj 
1996) to calculate model emission spectra, and of the GSI0 
and CFITSICEI external, auxiliary libraries. 

With respect to the structure of the code, the first unit 
is the most general, essential one, in which the emission as- 
sociated to the gaseous component in the simulation is gen- 
erated without including any additional specification, nei- 
ther about the observation to be simulated nor about the 
characteristics of the X-ray instrument, but the theoreti- 
cal emission model. Ideally, the spectral model assumed can 
be very complicated and account for more realistic descrip- 
tions of the emitting gas which are included in up-to-date 
sophisticated hydro-simulations, such as metallicity, chem- 
ical composition and sub-grid turbulent broadening of the 
emission lines. Conveniently, this first unit can be executed 
independently from the others just once per simulation out- 
put, enormously reducing the computational effort and time. 
The second unit of the code takes into account the geometry 
of the problem: the projection along a given line of sight is 
applied and photon energies are corrected for the Doppler 
Shift. The order of first and second unit basically reverses 
the approach of publicly available X-ray virtual telescopes 
and this has also the main advantage of delaying any limi- 
tation on spectral resolution as much as possible, preserving 
the possibility for any single photon to contribute to the 
final spectrum. 



2.1 Numerical hydro—simulation input 

PHOX requires as input the output of a hydrodynamical 
simulation. For the snapshots of the simulations we use, the 
redshifts and cosmological parameters are provided as well 
as the main quantities characterising the gas elements, such 
as density, temperature, metallicity, position and velocity. 
Currently, the simulator is best-suited to process simula- 
tions performed with SPH, particle-based codes, in partic- 
ular outputs of the N-body/SPH code GADGET-2 (e.g. 
ISpringel et all l200ll ; ISpringell l2005h . Nevertheless, the ap- 
proach is very general and the code can be easily adapted to 
process also grid-based simulations, for instance performed 
with Adaptive Mesh Refine ment code s like ART dKravtsovl 
1 19991 ; iKravtsov et~ail 120021 ). ENZ O (iNorman et ail l200l 
and references therein) or FLASH (jFrvxell et al.ll2000l ). In 
the description of the technique, we will preserve this gen- 
erality and refer to the gaseous component in terms of more 
general emitting volumes, namely "gas elements" , meaning 



The source code of the presented version of PHOX 
is made available by request to the authors (see 
http://www.mpa-garching.mpg.de/~kdolag/Phox/). We also 
plan to make the code available directly for anonymous download 
in the near future. 

2 See http://www.gnu.Org/s/gsl/ 

3 See http://heasarc.gsfc.nasa.gov/fitsio/ 



^ Hydrodynamical simulation output*^ 



UNIT1 



3_ 



For each gas element in the simulation 
output compute a model spectrum 



Calculate total number of photons from 
spectrum, assuming A fi4 , r ohsJid 



Populate the current spectrum with 
a distribution of photons 



Output files with photon packages. For each package 
position, velocity and a list of photon energies are stored. 



UNIT 2 




Choose l.o.s. direction 








Select subregion of interest 






[ Correct photon energies for Doppler Shift along l.o.s 1 






| Filter according to desired collecting area anc 


observing time | 









Photon list (arrival time, 
direction and energy) 



yes 



UNIT 3 

Convolve the simulated 
photon list with the 
response of a desired 
X-ray instrument 




Use any external 
software explicitely 
designed to simulate 
a specific X-ray 
instrument 
(e.g.,XISSIM) 



Synthetic X-ray event file 



Figure 1. Flow chart illustrating the main units of PHOX. The 
dashed box represents the possibility to couple PHOX with any 
external software specifically designed to simulate observations 
with a certain X-ray telescope. 



that these could be either particles, as in the application 
presented in the following Sections [3] and [4] or grid cells. 



2.2 Unit 1: generation of the box of photons 

The first, essential unit is devoted to convert the simulation 
snapshot from a box of gas elements into a box of photons. 
It is crucial to perform this step before the projection is 
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chosen and the emission is integrated along the line of sight 
in order to achieve and guarantee a very high spatial and 
spectral precision, accounting for the details of the three- 
dimensional structure of the source. 

Given the density, temperature and metallicity (or 
abundances) , the spectrum associated with each gas element 
of the simulation input is generated by means of the publicly 
available XSPEC package (v.l2jf|. In particular, since the ICM 
of galaxy clusters is our primary target, the code assumes ei- 



ther a MEKAL 



s our primary target, tnc code assumes ci- 
Mewe et all 1 19851; iKaastra fc Mewelll993l: 



iLiedahl et"allll995h or an APEC (|Smith et alj 120011 ) single 
temperature, thermal emission model for the hot plasma, 
with the possibility to add an absorption component give n 
by the WABS model (e.g. iMorrison fc McCammon] Il983h . 
Nonetheless, the code can be readily adapted to account 
for any desired emission model that is suitable for the con- 
sidered X-ray source (e.g. VMEKAL, VAPEC, BAPEC, BVAPEC, 
for more sophisticated descriptions of the ICM emission in- 
cluding also element abundances or velocity broadening for 
emission lines). At this stage, we do not convolve the model 
with any specific instrumental response but we use instead 
an identity matrix to provide very fine binning over a wide 
energy band. The choice of the model spectrum adopted, as 
well as of the energy limits and number of channels for the 
spectral binning, is left to the user, who can directly specify 
them in an input parameter file accepted by the code. Addi- 
tionally, a fiducial angular-diameter distance to the source 
(Da, fid) can be specified for convenience as an input param- 
eter, instead of calculating it accordingly to the cosmological 
redshift of the simulation output. 

The main calculation of this first unit consists of com- 
puting the emission associated to the gas component pro- 
cessing element per element, which even makes the code 
naturally suitable to parallelization. Our simulator directly 
populates the spectrum of each gas element with a distribu- 
tion of photons, without storing every spectrum or building 
a library of template spectra. Specifically, we calculate a 
cumulative spectrum and perform a Monte-Carlo-like sim- 
ulation to sample it with a total number of photons, N p h, 
determined by the total counts expected from the spectral 
model (in units of photons/s/cm 2 ) and the collecting area 
{A fid) and observing time (r bs,/id) specified by the user, so 
that 



N p h oc A/id x T obs ji d . 



(1) 



Ideally, the fiducial values for {A fid) and (r bs,/»d) have to 
be assumed in a convenient way such that the spectra are 
largely over-sampled. 

The photons are therefore characterised by energy (in 
the emitting element rest-frame), position and velocity. In 
order to reduce the amount of memory required to store the 
photon data, photons are organised in terms of packages, 
each package being associated to an emitting gas element. 
In this way, photon energies can be stored in a separate out- 
put file with respect to position and velocity, which in this 
first version of the code are equal for all the photons emitted 
by the same gas element and therefore stored just once per 
package. Moreover, the data are therefore naturally com- 
pressed, since positions and velocities need to be recorded 



See http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 



only for gas elements that indeed emit photons. The total 
number of photons produced from the whole gaseous com- 
ponent in the snapshot processed by PHOX is meant to be 
largely over-abundant in order to permit a dramatic reduc- 
tion in the following units because of geometrical selection, 
projection and instrumental response. 

The data generated from this basic part of PHOX rep- 
resent a "cube of photons" associated to the input hydro- 
simulation. From this stage on, the simulator has all the 
information required to proceed with the synthetic observa- 
tion and the original input is no longer needed. Therefore, 
it is worth to remark that this first unit can be processed 
once per output independently of the specific study to be 
performed. 



2.3 Unit 2: projection and preparation for the 
observation 

PHOX second unit takes into account the geometry of the 
observation to be simulated. It accepts the output data 
produced by the first unit and a parameter file with user- 
specified parameters for the sub-region of the cube that has 
to be selected, direction corresponding to the line of sight 
(l.o.s.), real collecting area and time for the mock observa- 
tion. 

After selecting geometrically the photons, the energies 
have to be corrected for the Doppler Shift in order to ac- 
count for the line-of-sight velocity of the gas element that 
originally emitted them. In particular, the correction from 
the emitted to the observed frame is expressed as 



e: 



ph 



E 



p h 



1-/3 
1 + /3' 



(2) 



where /3 = vi. . s ./c and Vi. . s . is the velocity component 
along the line of sight. For the sake of simplicity, the line 
of sight is assumed to be aligned with the z-axis. Neverthe- 
less, any different, desired direction only requires a simple 
rotation of the package positions and velocities before run- 
ning this second unit of the code. 

The sample of photons obtained has also to be further 
filtered according to the specified collecting area, A, and 
a realistic exposure time, r b s , chosen to simulate the ob- 
servation. Accordingly to Eq.Q] this is done by calculating 
the re-scaling factor between these values and the fiducial 
quantities assumed in the first unit as 



fak 



A T bs 



Afid Tohs.fid 



r> 2 

u A,fid 

n 2 



(3) 



and by assuming it as a probability factor for each photon to 
be actually observed. Eq.[3] explicitly accounts for the possi- 
bility of re-sacling the number of observable photons by the 
angular-diameter distance to the source, Da, whenever this 
parameter differs from the fiducial value, Da, fid, adopted in 
Unit 1. 

Lastly, the photon list has to be stored in the most 
convenient format either to be convolved with the techni- 
cal response of the specific instrument or to interface with 
an external software dedicated to simulate specific X-ray 
satellites. 

This second unit already takes into account the charac- 
teristics of the study to be performed but no longer requires 
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the original simulation data, processed by Unit 1. Therefore, 
any change in the specification of the mock observation can 
be easily included in this post-processing phase. As such, a 
straightforward consequence is the unique possibility to in- 
vestigate the same astrophysical source from many different 
line of sights, handily processing with Unit 2 the cube of 
virtual photons obtained with the first unit. 

2.4 Unit 3: simulating the observation 

At the third and last stage the mock X-ray observation is 
completed by eventually considering a real telescope. The 
photon list obtained from the second unit is convolved with 
the technical characteristics of a specific instrument, defined 
by the redistribution matrix file (RMF) and the ancillary 
response file (ARF). Such process produces final event files 
which satisfy the standards of real X-ray observations, so 
that they can be analysed with the same procedures and 
tools used for real data. 

In the most general perspective this last unit assumes 
the RMF and ARF files supplied by the user and it is not 
constructed a priori for a specific instrument. Particular at- 
tention, in general, has to be payed to the normalization of 
the effective area defined by the ARF file of the instrument 
with respect to the effective area assumed during the pro- 
jection phase, in order to avoid unphysical overabundance 
of observed photons. 

The third unit of the simulator, being independent 
of the others, can be conveniently replaced by any de- 
sired X-ray instrument simulator, like the sophisticated tool 
XlSSIlJf], which has been developed especially to obtain syn- 
theti c observations of th e Suzaku X-ray Imaging Spectrom- 
eter (|lshisaki et al.ll2007h . 



3 SYNTHETIC OBSERVATION OF THE 
COSMIC WEB: TWO CASES OF STUDY 

We present here an optimal application of PHOX to the 
hydro-numerical simulation of a filament-like structure. The 
large-scale region provides indeed a study case to test the re- 
markable gain in computational cost and time for the post- 
processing of large data sets; also, the cluster-like haloes re- 
siding in this simulated filament offer various science cases 
to perform the spectral analysis of the mock X-ray observa- 
tions. 



3.1 The simulated region 

The high-resolution simulated region containing the 
supercluster-like structure extends for about 50 x 50 x 
70Mpc 3 and consists of 27 haloes at present epoch, four 
of them being massive cluster-size haloes. This filamentary 
region has been o riginally extracted from a cosmological, N- 
body simulation (| Jenkins et al.ll200ll ; lYoshida et aTl200 ll ) of 
a ACDM universe with fi = 0.3, Q b = 0.044, a 8 = 0.9 
and h = 0.7, within a box of 479/i -1 (comoving) Mpc a 
side. Using the "zoomed initial condition" (ZIC) technique 



Table 1. Typical computational time scales and output storage 
memory for the hydrodynamical numerical simulation and differ- 
ent units of PHOX. As a representative case, the data reported 
here refer to the science application t presented in Section [3] 



Sim 




PHOX 






Unit 1 


Unit 2 


Unit 3/XISSIM 




~lh 40min* 


~few sec 


~few min 


1.8 Gb 


827 Mb 


176 Mb (CL2) 
60 Mb (CL1) 


~100 Mb (CL2) 
~35 Mb (CL1) 



* The run reported here has been performed by binning in temperature 
the particles (assuming A(kT) ~ beV) and calculating a spectrum 
per bin, rescaling afterwards the expected total number of photons by 
the specific normalization calculated for each particle in the temper- 
ature bin. Since no varying mctallicitics have been considered in this 
run, this approach has been used to further reduce the computational effort. 

f The results arc obtained for serial runs performed on standard 
work station (2300 MHz, AMD Opteron). 

(|Tormen et all Il997l) . the overdense region has been re- 
simulated at higher resolution and a hydrodynamical run 
with these new initial conditions has been performed with 
the TreePM/S PH code GADGET-2 l|Springel et all l200ll ; 
Springcl 2005), including cooling, star formation and feed- 
back from supernova winds. The simulation output contains 
approximately 1.2 x 10 7 gas particles and 1.7 x 10 7 DM 
particles, at z = 0, and the final resolution is tudm = 
1.30 x 10 9 /i _1 M Q and m gas = 1.69 x 1O 8 /i _1 M , for the 
DM and gas particles, respectively. For an extensive and de- 
tailed description of the si mulation run and an alysis we refer 
the reader to the paper bv lDolag etail (|2006r i. 

3.2 Suzaku mock observations with PHOX 

We applied PHOX to the whole simulation output corre- 
sponding to redshift z = 0.07 and obtained from the first 
unit a list of all the photons emitted by the sufficiently-hot 
gas particles in the simulated region. In order to compute the 
model spectra, we adopted an APEC model with fixed metal- 
licity, set for simplicity to Z = O.3Z0 for all the particles, 
and a ssumed the solar abundances by lAnders fc Grevessd 
(1989). The value assumed for the Galactic equivalent col- 
umn density for the WABS absorption model was Nh 
7 x 10 20 cm -2 . The values for fiducial collecting area and 
exposure time were initially set to 2000cm 2 and lMs, re- 
spectively. This first stage was the most computationally 
demanding (see Table l3T2"l for time scales and memory re- 
quirements referring to the application presented here) and 
permitted to generate approximately 1.8 x 10 8 photons, for 
the whole simulation output. 

In Unit 2, we assumed to observe the large-scale re- 
gion from a physical distanc^fl of 313.9Mpc, which corre- 
sponds to the luminosity distance at the snapshot redshift 
z = 0.07. Similarly, with the given cosmological parameters, 
the angular-diameter distance at this redshift is 274.4Mpc, 
which gives a physical scale of 1.33kpc/arcsec. In this pro- 
jection phase, we select two specific sub-regions of the sim- 



6 We set the observer on the positive part of the z— axis, assumed 
5 See http://heasarc.nasa. gov/docs/suzaku/prop_tools/xissim/xissim_utogscb±h)^iied with the l.o.s.. 
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3Mpc 




Figure 2. Filament-like region and zoom onto the haloes CL1 and CL2, at z = 0.07. Left-hand panel: projection along the 2-axis of 
the emission— weighted temperature, in a logarithmic color scale. The map is 50Mpc a side and encloses the high-resolution region of the 
box containing the filament. The 3 X 3Mpc 2 regions centered on CL1 and CL2, respectively, are marked with the green squares. Central 
column: zoom onto CL1 (top panel) and CL2 (bottom panel). The maps are 3Mpc a side and show the logarithm of the emission-weighted 
temperature, projected along the 2-axis. The -R500 region (cyan circle) and the XIS FOV (magenta dashed square) are overplotted for 
both haloes. Right-hand column: theoretical spectra obtained from Unit 1 and 2 of PHOX for the R500 region of both haloes. 



ulation box, containing two cluster-like haloes among the 
four most massive ones. In particular, we consider for the 
present analysis a regular, massive cluster (hereafter, CL1) 
with A/500 = 3.97 x 1O 14 M and #500 = 1070.5kpc and 
a massive, disturbed system with M500 = 9-24 x 1O 14 M0 
and f?5oo = 1417. 3kpc (hereafter, CL2). The region geomet- 
rically selected around both haloes is a cylinder along the 
z-axis, enclosing -R500 in the xy plane. Planning to simu- 
late a Suzaku observation, the collecting area assumed for 
the mock observation equals the physical area of the Suzaku 
X-ray Imaging Spectrometer (XIS), A = 1152.41cm 2 , and 
the exposure time is T ts = 500ks, in order to have good 
statistics. This allows to extract 3.1 x 10 6 photons for CL1 
and 9.2 x 10 6 photons for CL2. 

In the left-hand panel of Fig.[2] we show the filament- 
like structure at z = 0.07, projected along the 2-axis (i.e. 
the observation l.o.s.). The map represents the emission- 
weighted temperature of the simulated filament enclosed in 
the high-resolution region. Overplotted in green, we mark 
the zoom onto the selected haloes, which are shown in the 
central column of Fig. [2] The two panels contain in fact the 
3 x 3Mpc 2 maps of the emission-weighted temperature pro- 
jected along the l.o.s., for CL1 and CL2. For comparison, 
we overplot also the -R500 region (cyan circle) and the XIS 
field of view (magenta dashed square). Instead, in the right- 
hand column of Fig.[2] we show the ideal integrated spectra 
for CL1 and CL2, as obtained from Unit 1 and 2 after se- 
lecting the photons coming from the region enclosed within 
-Rsoo- These spectra are obtained purely through a regular 
binning of the photon energies with bins of AE = O.OOlkeV, 
without including any sensitivity or effective area. 



For the purpose of this paper, we present PHOX by tak- 
ing a dvantage of the public package XISSIM (Ishisa ki et al.l 
2007), designed to simulate observations with the XIS spec- 
trometer on board the Suzaku X-ray satellite. By adapting 
the second unit of our code to produce the output photon 
list in a format best-suited for XISSIM, we test our novel 
technique and demonstrate its capabilities through a simple 
science application. 

Specifically, we use XISSIM to process the simulated 
photons generated by PHOX with the real responses and 
calibration characteristics of different XIS CCDs. Given 
the high statistics offered by our simulation, we decide to 
consider the back-illuminated (BI: XIS-1) sensor and the 
two combined, front-illuminated CCDs (FI: XIS-0,-3) sepa- 
rately. The combination of two 500ks observations with the 
XIS-0 and XIS-3 detectors additionally improves the statis- 
tics for the final FI spectrum. 

The event files obtained for the two cluster-like haloes 
of the filament, CL1 and CL2, without any addition of phys- 
ical background emission, were then analysed following the 
standards of X-ray data analysis. 



3.3 Spectral Analysis 

In Fig.0 we show the simulated photon images (left-hand 
column) of the 500ks XIS observations, with the BI and FI 
sensor, for CL1 and CL2. The images, as well as the spectra, 
were extracted from the original event files using Xselect 
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(v.2.4) from the FTOOL^I package (|Blackburnlll995h . and 
correspond to the central region in both clusters. In partic- 
ular, the Suzaku XIS covers a field of view (FOV) of roughly 
18' x 18', which corresponds to a physical scale of 1436. 6kpc 
at the distance of our sources and encloses therefore the re- 
gion within ~ 0.67i?5oo of CL1 and ~ 0.51-Rsoo of CL2 (we 
refer to the central panels of Fig.[5]for a visual representation 
of the R500 region and the XIS FOV in the two haloes). 

Simulations of cluster-like haloes provide the possibil- 
ity of precisely knowing the intrinsic dynamical and thermal 
properties of the ICM. From a visual inspection of the tem- 
perature maps shown in Fig. [2] (central column), the differ- 
ent structure of the two haloes studied is already evident. 
As we will discuss in more detail in Section [4] (see Fig.Q, 
the further investigation of the emission measure (hereafter, 
E.M.) distribution as function of temperature, for the gas 
particles in the simulation that reside in the FOV of Suzaku 
both in CL1 and CL2, let us unveil the halo intrinsic ther- 
mal structures. Especially, CL1 is mainly dominated by one 
temperature component, although the E.M. distribution is 
not very narrow, whereas the second halo, CL2, clearly has a 
complex thermal structure, which cannot be well described 
by a single temperature. 

Therefore, we fit a multi-temperature model to the BI 
and FI XIS spectra, restricting our analysis to the 0.2 — 
10 keV energy band and requiring a minimum of 5C0 counts 
per en ergy bin. The spectral fit was performed with XSPEC 
12.6.0 |Arnaudlll996t ). The response files (RMFs) used in the 
fit were the same as in the run of XISSIM, while we generated 
the ARF files for ea ch halo and detect or by means of the 
ftool XISSIMARFGEN jlshisaki et al.ll2007t) . 

For our spectral modelling of the simulated data, we 
considered a WABS x (APECi + APEC 2 + APEC 3 + APEC 4 + APEC 5 ) 
model, describing each component by an APEC model, as 
assumed in Unit 1 for the emission associated to each 
gas element. The equivalent hydrogen column density for 
the Galactic absorption model WABS was frozen to Nh = 
7 x 10 20 cm -2 , as in the original run of Unit 1. For anal- 
ogous reasons, we adopted fixed values for redshift, z = 
0.07, and metallici t y Z — 0.3ZQ, always assuming the 



lAnders fc Grevessd (l989T ) solar abundances 

In both haloes, we assigned to the first temperature 
component, T*, a value fairly close to the emission-weighted 
temperature estimated from the gas particles in the simula- 
tio n. Then, followi n g a st rate gy similar to w h at is s uggested 
by iPeterson et alj 120031 ) or IKaastra et all l |2004 ). we as- 
signed to the other components temperatures of 2T* and 
|T«, jT*, gT*, for the cooler ones. The five normalizations 
of the APEC components where free parameters in the fit. 

The spectra and best-fit models are shown in the right- 
hand column of Fig. [3] for the two haloes and the two XIS 
sensor kinds. The results of the spectral fit are presented 
in Table EOl The five temperatures, frozen in the fit, are 
different for CL1 and CL2 and therefore reported in the table 
for clarity reasons. Overall, the resulting Xred i s always very 
good, meaning that the five-temperature model provides a 



fair description of the thermal structure for both clusters in 
the central observed region. 



4 RESULTS: RECOVERING THE EMISSION 
MEASURE DISTRIBUTION 

The interesting purpose of our multi-temperature modelling 
is the reconstruction of the thermal structure of the emitting 
ICM. This is pursued by directly comparing the distribution 
of the best-fit normalizations at the corresponding, fixed 
temperatures with the expected values calculated from the 
simulation. Dealing with simulated data, we can faithfully 
test the results of the spectral analysis against the underly- 
ing solution in a straightforward way, which is not possible 
for real data. 

By definition, the normalization of each APEC compo- 
nent is 



K = 10~ 14 E.M./[47rDi(l + z) 2 ] cm" 



(4) 



where Da is the angular-diameter distance to the source 
and z is the redshift. This implies that the normalization is 
directly proportional to the E.M. of the hot emitting plasma, 
defined as 



E.M. 



n e riHdV, 



(5) 



7 See http://heasarc.gsfc.nasa.gov/ftools/ 

8 The spectra were rebinned using the grppha routine from the 
FTOOLS package. 



where n e and uh are the electron and hydrogen number den- 
sities of the plasma, respectively. In order to map the distri- 
bution of E.M. in the simulated haloes to the corresponding 
K values at fixed temperatures, the theoretical temperature 
binning has to be determined in the most faithful way pos- 
sible. 

The shape of the E.M. distribution extracted from the 
simulation is presented for the two haloes in the left-hand 
panels of Fig. [4] This has been calculated considering all 
the gas particles of the simulated haloes, for which the pro- 
jected position on the xy plane is enclosed in the Suzaku 
field of view. A lower temperature cut has been artificially 
set to Tmin = 0.2keV, which is the lower boundary of the 
energy band used in the spectral fit. Nevertheless, this con- 
tinuous distribution has to be reduced to the five fixed 
temperatures adopted for the spectral fit. The most con- 
sistent and independent way of choosing the boundaries 
between the different temperature components, is by re- 
quiring that the value of the E.M.-weighted temperature 
for each bin equals the fixed value adopted for the fitting 
(w e follow here an appro ach analogous to the one adopted 
by IKaastra et al.1 ([2004)). Specifically, we iteratively ad- 
just the division boundaries until the match between the 
E.M.-weigthed temperature in the bin and the fixed value 
is reached, starting from the lower-temperature component 
and proceeding towards the higher ones. Finally, the total 
E.M. calculated in each bin can be converted into an ex- 
pected value for the APEC normalization via Eq.|?] 

The direct comparison between simulation and mock- 
observation results is presented in the right-hand panels of 
Fig.FJ] for both CL1 and CL2. The original E.M. distribution 
for the two haloes is plotted in the left-hand panels of Fig. 31 
where it is clear that none of them is strictly isothermal, al- 
though the CL1 cluster has a narrower thermal distribution 
with respect to CL2. The theoretical expectations for the 
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1000 2000 4000 Energy (keV) 



Figure 3. Simulated Suzaku images (left-hand column) of CL1 and CL2 and corresponding spectra (right-hand column). Also the 
models fitted to the spectra are shown. For each halo the BI XIS and the combined FI XIS sensors are presented. 
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Table 2. Multi-temperature best— fit results for Suzaku XIS simulated spectra of haloes CL1 and CL2. The model adopted is WABS X 
(APECi + APEC2 + APEC3 + APEC4 + APEC5), where only the normalizations Ki are free in the fit. 



XIS 


K* 


K 2 


K 3 


K 4 


K s 


X 2 /d.o.f. 


CL1 




(fcTi = 0.69keV) 


(kT 2 = 1.38keV) 


(kT 3 = 2.75keV) 


(kT 4 = 5.5keV) 


(fcT 5 = llkeV) 




BI 
FI 


0.21(±0.04) x 10~ 3 
0.27(±0.04) x lCr 3 


0.15(±0.16) x 10" 3 
0.12(±0.14) x 10~ 3 


2.03(±0.60) x 10~ 3 
2.00(±0.44) x 1CT 3 


8.17(±0.63) x 10~ 3 
8.12(±0.43) x 1CT 3 


1.38(±0.54) x 10" 3 
1.47(±0.36) x lCT 3 


1693.2/1700 
1713.2/1998 



CL2 



(fcTi = IkeV) 


(kT 2 = 2keV) 


(fcT 3 = 4keV) 


(fcT 4 = 8keV) 


{kT 5 = 16keV) 




BI 0.99(±0.09) x 10" 3 
FI 1.16(±0.09) x 10- 3 


0.0(±0.91) x 10- 3 
0.0(±0.70) x 10- 3 


7.87(±1.70) x 10~ 3 
7.49(±1.21) x lO" 3 


16.44(±2.28) x 10~ 3 
16.14(±1.52) x lO" 3 


7.87(±1.41) x 10" 3 
8.32(±0.90) x 10~ 3 


2762.0/2167 
2557.3/2435 



*The normalization of an APEC component is defined as K = 10 14 / rienndV ~/[i-KD 2 A (l + z) 2 ] cm 5 , where Da is the 
angular-diameter distance to the source. 
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Figure 4. Left-hand column: emission measure (E.M.) distribution as a function of temperature for gas particles contained in the FOV 
of Suzaku XIS (projected along the z-axis), for haloes CL1 and CL2 at redshift z = 0.07. Right— hand column: comparison between 
best-fit values of APEC (triangles) and BAPEC (—1, diamonds and -2, squares) normalizations and corresponding theoretical values 
(asterisks), computed from the E.M. distribution of emitting gas particles. Different colors for the best-fit values refer to spectra obtained 
with different XIS detectors: BI XIS (red), and FI XIS (sum of XIS-0 and XIS-3, green). In both columns, the upper row refers to CL1 
and the lower row refers to CL2. 
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fixed temperature components are plotted as black aster- 
isks on the right-hand panels, where the best-fit values are 
reported as well. 

Interestingly, the agreement between simulated and 
"observed" values is quite good and the overall shape of 
the thermal distribution in both clusters is recovered rea- 
sonably well. This result is found to be true for both the 
BI and FI XIS sensors. Moreover, the differences between 
the thermal structures of the two haloes are reflected by a 
different width and extension of the E.M. distribution over 
temperature, which is definitely broader in the CL2 case. 

We notice that the dominant component at kT^ — 
5.5keV is underestimated by the fit for the halo CL1. 
Nonetheless, it is worth to notice that the determination of 
the most consistent temperature binning of the theoretical 
E.M. distribution is expected to be sensitive to the defini- 
tion of temperature assumed to match the temperature of 
the thermal component in the spectral fit. For instance, it 
is known that there is a discrepancy between the emission- 
weighted temperature, as usually calculated in the analysis 
of numerical simulations, and the projected spectroscopic 
temper ature obtained from X -ray observations (see, for in- 
stance, iMazzotta et al.l [2004) . In this respect, our choice 
for the theoretical temperature calculated from the simu- 
lation, which is weighted by the E.M., is believed to be the 
most consistent. However, we also tested the results against 
a different definition of temperature, namely the emission- 
weighteojf] one, which resulted in a slightly different set of 
theoretical values for the five APEC normalizations. Even in 
this case, our conclusions are not affected in any significant 
way and the overall shape of the underlying E.M. distribu- 
tion is equally well traced by the multi-component spectral 
fit. 

In the right-hand plots of Fig. [3] we also show two al- 
ternative five-temperature models adopted to fit the syn- 
thetic spectra. Specifically, we perform two spectral fits 
with a multi-temperature model consisting of five BAPEC 
components. The BAPEC spectral model is a velocity- and 
thermally-broadened emission spectrum, where the velocity 
broadening is parametrized by a gaussian velocity disper- 
sion a v . In fact, the technique implemented in our simulator 
allows to preserve the very high energy resolution of the 
integrated spectrum till the final convolution with the in- 
strument. Since every photon energy carries the imprint of 
the emitting particle velocity, the spectral features like metal 
emission lines could include a non-thermal broadening due 
to peculiar non-thermal motions of the gas particles of the 
ICM and the velocity-broadened APEC is a valid choice to 
take this into account. In the two cases explored with the 
WABS X (BAPECi + BAPEC 2 + BAPEC 3 + BAPEC 4 + BAPEC 5 ) model, 
every frozen parameter was assumed in the same manner as 
in Section 13.31 and the five BAPEC normalizations were let 
be free. In addition, in the first model (BAPEC-f), we fit for 
another free parameter, which is the velocity dispersion a v 

9 The emission-weighted temperature is defined as 
_ fn*A(T)TdV 
Jn 2 A(T)dV ' 

where n is the gas density and A(T) is the cooling function. As 
usually done, we consider here A(T) oc T3, assuming the gas to 
emit mainly via Bremsstrahlung. 



linked across the five components, whereas in the second 
case (BAPEC-2) both a v and the redshift are left free in the 
fit. As shown in the plots, the results for the normalizations 
are definitely consistent with the standard case, where the 
five thermal components are described by the APEC model. 
However, the best-fit values for er„, in particular, turn out to 
have very large errors, meaning that this additional param- 
eter is not reliably fitted. We interpret this as an evidence 
that the energy resolution of the Suzaku XIS spectrometer 
is not sufficiently high to capture the details of the emission 
lines and therefore a velocity-broadened model like BAPEC is 
not improving the fit to the synthetic data. 

Moreover, the results look generally stable against 
metallicity, which could in principle compete with normal- 
ization in order to reproduce the spectral features associated 
to metal emission lines. In particular, we tested the fit with 
five APEC models, both for CLI and CL2, allowing for the 
metallicity value to be free, although linked across the five 
components, in addition to the five normalizations. Gener- 
ally, this converges to values within few percents from the 
expected metallicity of Z = 0.3Zq, used in the simulation 
of the X-ray observation. Comparing to the results of the 
standard fit reported in Table [<01 the APEC best-fit normal- 
izations, for initial metallicity values set to [0.2, 0.3, 0.5]Zq , 
are found to be consistent with the original values, within 
the error bars associated to them. 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a novel X-ray photon simu- 
lator, named PHOXpl. dedicated to obtain synthetic X-ray 
observations from the output of hydrodynamical numerical 
simulations. 

The development of PHOX has been strongly motivated 
by the ultimate aim of comparing the output of hydrody- 
namical simulations of cluster-like objects to real X-ray ob- 
servations in the most faithful and reliable way possible. In 
fact, X-ray observations still provide us with increasingly 
detailed data that require a deep understanding of the un- 
derlying structure of galaxy clusters. Numerical simulations, 
on the other side, continuously improve in describing not 
only the dark matter but also the complex physical processes 
governing the baryonic component. Therefore, a strategy de- 
voted to combine the advances and results in both observa- 
tions and simulations can be promisingly successful. 

In the past, software packages like X-MAS 
|Gardini et all |2004 ) or XIM (|Heinz fc Briiggenl l2009h 
have been developed to produce synthetic observations 
of astrophysical sources simulated with particle- and 
grid-based codes, respectively. The existence of such tools 
is definitely complementary to our simulator in order to 
test the robustness of comparison studies between state- 
of-the-art hydrodynamical simulations including different 
physical phenomena and real observations. With respect to 
other codes, the novelty of our virtual telescope resides in 

10 The source code of the presented version of PHOX 
is made available by request to the authors (see 
http://www.mpa-garching.mpg.de/~kdolag/Phox/). We also 
plan to make the code available directly for anonymous download 
in the near future. 
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the original method implemented. In particular, the order 
of the three main units that constitute the code offers a 
significant gain in computational cost and time, making 
PHOX extremely flexible and useful to process hydro- 
numerical simulations of several astrophysical objects, at 
various spatial scales, as well as simulated light-cones. 
Moreover, this method uniquely guarantees extremely high 
spatial and energy resolution. The more demanding phase is 
indeed the first one, where the output of a hydrodynamical 
simulation is processed and the X-ray emission of each gas 
element (either particle or grid cell) is simulated by using 
the spectral models available in the X-ray software XSPEC. 
Essentially, we perform the sampling of each spectrum with 
an expected number of photons, and we immediately store 
the sample of photons insted of constructing a library of 
template spectra to be populated afterwards. As such, the 
spectral model chosen to mock up the emission from the 
gas elements, regardless of its complexity, can be explicitly 
specified and suited to represent any desired astrophysical 
source. Also, the relatively reasonable requirements in terms 
of computational time and memory, make out of PHOX 
an ideal tool to process forthcoming large cosmological, 
hydrodynamical simulations as well as the usual smaller 
ones. 

Concerning galaxy clusters, PHOX currently uses two 
possible models for the X-ray emission of the ICM (namely, 
MEKAL and APEC) and accounts for the particular density 
and temperature of the simulated gas elements. However, 
the technique implemented offers the chance to account also 
for different metallicities, or even metal compositions, across 
the emitting gas elements, without increasing the memory 
requirements significantly. 

Furthermore, PHOX postpones the projection along the 
line of sight and the convolution with a real instrument re- 
sponse till after the generation of the discrete sample of 
photons, permitting to each of them, in principle, to con- 
tribute to the final integrated spectrum. Requiring a high 
energy resolution in the model spectrum associated to every 
gas element, we can ideally combine the potential of PHOX 
with the characteristics of upcoming X-ray satellites, like 
1X0/ ATHENA or ASTRO-H, in order to predict the de- 
tectability of fundamental spectral features. This could be 
of great help, for instance, in constraining the ICM non- 
thermal motions (e.g. turbulent, bulk, rotational motions) 
via the velocity broadening of heavy-ion emission lines. In 
general, the theoretical modelling of X-ray-emitting sources 
is a powerful way to predict and optimize the capabilities of 
future X-ray missions and to help in the interpretation of 
the data when the instruments are in operation. 

As an interesting application of this new virtual 
telescope, we focused on the hydro-numerical simulation 
of a large-scale region, originally performed with the 
TreePM/SPH code GADGET-2 and then re-simulated in- 
cluding cooling, star formation and feedback from super- 
nova winds. This filament -like structure is particularly wor- 
thy for presenting the simulator in each part and espe- 
cially for demonstrating the advantages of our approach. 
From a purely technical perspective, the greatest gain was 
achieved in terms of computational cost since the entire, 
large data cube contained in the simulation output was 
processed with PHOX Unit 1 only once, providing an out- 
put of relatively compacted photon data (see Table l3T2"j) . In 



fact, this first phase permitted to obtain the X-ray pho- 
tons virtually emitted by all the gas particles in the fila- 
ment. At a later stage, we restricted our analysis to two 
cluster-like haloes with different thermodynamical struc- 
ture, CLI and CL2, which were geometrically selected from 
two different regions of the simulation box. Conveniently, 
the projection and observation units of PHOX could be 
run separately for each halo processing the same output 
generated with Unit I and we finally obtained synthetic 
Suzaku images and spectra (Fig. [3]). The convolution with 
the response of the X-ray Imaging Spectrometer (XIS) on 
board Suzaku, via the public tool XISSIM, provided data 
that could be analysed with the standard procedures of 
X-ray analysis. An important aim of this application was 
the reconstruction of the intrinsic emission measure distri- 
bution by the multi-temperature fitting of the simulated 
spectra. The goal of setting spectroscopic constraints on 
the multi-phase ICM of galaxy clust ers was pursued in X - 
ray observational stu dies reported by Peterson et ail {2003), 
iKaastra et all l|2004 ) and lFabian et all (120061 1. although for 
real clusters it is intrinsically impossible to know the un- 
derlying multi-phase structure unambiguously. Instead, nu- 
merical simulations uniquely allow to model and predict the 
theor etical temperature distributio n of the ICM (e.g., stud- 
ied in Ka wahara et al.ll2007l . |200ST ) and constitute therefore 
the ideal possibility to test the power of the method used to 
reconstruct the structure of the mock data. The direct com- 
parison between the X-ray data and the theoretical thermal 
distribution of the input simulation is found to be overall 
faithful and robust for both CLI and CL2 (see Fig.QJ. More- 
over, the different thermal structure that distinguishes the 
two selected haloes is remarkably highlighted by the multi- 
temperature fitting. This result promisingly suggests that 
the temperature distribution of the ICM can be in princi- 
ple traced and unveiled by spectral analysis of X-ray ob- 
servations. From the additional investigation of velocity- 
broadened fitting models (like the BAPEC model) for the 
Suzaku spectra, we can also conclude that a higher spec- 
tral resolution is required in order to reliably constrain the 
ICM velocity structure. In general, it is important to notice 
that the application discussed here does not include other 
emission components (e.g. a physical X-ray background) in 
addition to the ICM one and represents therefore an ideal 
case study. 

With regard to technique, the main feature distinguish- 
ing PHOX from other virtual X-ray telescopes, such as 
XMAS and XIM, is the explicit creation and storage of a 
3-dimensional box of photons from the raw simulation data 
cube itself. This step is performed before and without any 
reference to the observation process and its result is stored 
as a separate data product, which represents the ideal X-ray 
emission from all the gas volumes in the simulated region, 
allowing us to employ Unit 1 of PHOX to create a photon 
database associated to the simulation. In Unit 2 this photon 
cube is projected onto a particular observational plane and 
the photons are sampled according to the chosen aperture 
and observation time, as described in Eq.(3] With the photon 
box pre-calculated, this step is computationally trivial and 
can be performed very efficiently. As such, it could in fact be 
made available through an on line web service. The result of 
this step is a photon list and is processed in this paper by 
means of the Suzaku simulator XISSIM. Similarly, we believe 
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it should be possible for most such virtual telescopes to be 
adapted to the same format without undue difficulties. 

This leads us to an interesting modular design, with a 
photon simulator producing a list of photons directly from a 
hydrodynamical simulation, and a virtual telescope receiv- 
ing this list, and convolving it with the characteristics of a 
real X-ray instrument to produce a photon event file that 
can be directly compared in all its aspects to the analogous 
files resulting from real observations. This design presents 
a proper separation of concerns between theorists and ob- 
servers, where each module is clearly independent from the 
others and can therefore be implemented and used sepa- 
rately. In a larger perspective, the well specified model for 
the data shared and exchanged between these units could 
represent a first step towards a standardised, interoperable 
way to publish simulation results in a way that is suitable 
for investigation by observers. In particular, this is very suit- 
able for the kind of approach pursued by the International 
Virtual Observatory Alliance (IVOA0), and we indeed plan 
to deploy a web service that provides access to the results 
discussed in this paper in that contexlF 2 !. 
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